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Abstract 

Nonlinear non-stationary equation describing evolution of weakly curved pre- 
mixed flames with arbitrary gas expansion, subject to the Landau-Darrieus 
instability, is derived. The new equation respects all the conservation laws 
to be satisfied across the flame front, as well as correctly takes into account 
influence of vorticity, generated in the flame, on the flame front structure and 
flame velocity. Analytical solutions of the derived equation are found. 
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I. INTRODUCTION 



Description of premixed flame propagation is, in essential, the description of develop- 
ment of the Landau-Darrieus (LD) instability of zero-thickness flames. Given arbitrary 
flame front configuration, its evolution is determined by the exponential growth of unstable 
modes, eventually stabilized by the nonlinear mode interaction. By themselves, the nonlin- 
ear effects are not sufficient to stabilize flame propagation, since the spectrum of unstable 
perturbations of a zero-thickness flame is unbounded. In many cases, however, an upper 
bound for the mode wavenumber is provided by the heat conduction - species diffusion 
processes in the flame, which govern the evolution of short- wavelength perturbations 0. 
The flame propagation is thus described as the nonlinear propagation of interacting modes 
of zero-thickness front, with an effective short- wavelength cut-off described at the hydro- 
dynamic scale by means of an appropriate modification of the evolution equation and the 
conservation laws at the flame front [Q. 

In this purely hydrodynamic formulation, the flame dynamics is governed essentially by 
the only parameter - the gas expansion coefficient 9, defined as the ratio of the fuel density 
and the density of burnt matter. Unfortunately, in general it is very difficult to reduce the 
complete system of hydrodynamic equations governing flame dynamics to a single equation 
for the flame front position. This is mainly because the nonlinearity of flame dynamics 
cannot be considered perturbatively. For instance, it can be shown that in the regime of 
steady flame propagation, the flame front slope can be considered small only if the gas 
expansion is small [6 —>■ 1), while for flames with 6* = 6 -^ 8 it is of the order 2-^3 (discussion 
of this issue can be found in Ref. [H). At the early stages of development of the LD- 
instability, however, the perturbation analysis is fully justified, and the equation describing 
nonlinear propagation of the flame front can be obtained in a closed form. Within accuracy 
of the second order in the flame front slope such an equation was obtained by Zhdanov and 
Trubnikov (ZT) ^j, without taking into account the influence of the effects related to finite 
flame thickness, mentioned above. The latter were included in the ZT-equation ad hoc by 
Joulin [^]. 

Concerning the ZT-equation and its modification, we would like to note the following. 

1) Although this equation respects all the conservation laws to be satisfied across the 
flame front, it was derived on the basis of a certain model assumption concerning the flow 
structure downstream. Namely, it was assumed that the velocity field can be represented 
as a superposition of a potential mode and an "entropy wave", so that the pressure field 
is expressed through the former by the usual Bernoulli equation. The generally nonlocal 
relation between the pressure and velocity fields is thus rendered local algebraic, which 
allows simple reduction of the system of hydrodynamic equations to the single equation for 
the flame front position. Being valid at the linear stage of development of the LD-instability, 
this model assumption is, of course, unjustified in general. 

2) It was assumed in the course of derivation of the ZT-equation that not only the front 
slope, but also the value of the front position itself is a small quantity. As a result, the ZT- 
equation turns out to be non-invariant with respect to space translations in the direction of 
flame propagation. Practical consequence of this assumption is the unnecessary limitation of 
the range of validity of the equation: the space and time intervals should be taken sufficiently 
small to ensure that the deviation of the front position from the initial unperturbed plane 
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is small everywhere. 

The purpose of this paper is to show that the only assumption of smallness of the flame 
front slope is actually sufficient to derive an equation describing the nonlinear development 
of the LD-instability to the leading (second) order of nonlinearity. Surprisingly, this equation 
turns out to be of a more simple structure than that of the ZT-equation. This simplification 
is due to existence of a representation of the flow equations at the second order, which can 
be called transverse. In this representation, the system of equations governing the flame 
propagation can be brought into the form in which dependence of all dynamical quantities 
on the coordinate in the direction of flame propagation is rendered purely parametric. 

Let us consider the question of existence of this representation, and more generally, 
the meaning of the weak nonlinearity expansion, in more detail. First of all, the following 
important aspect of the problem should be emphasized. The curved flame propagation is an 
essentially nonlocal process, in that the presence of vorticity produced in the flame implies 
that the relations between flow variables downstream generally cannot be put into the form 
in which the value of one variable at the front surface can be expressed entirely in terms of 
other variables taken at the same surface. For instance, the value of the pressure field at the 
flame front depends not only on the gas velocity distribution along the front, but also on its 
distribution in the bulk. This means in turn that the equation describing the flame front 
evolution cannot generally be written in a closed form, i.e., as an equation which expresses 
the time derivative of the front position via its spatial gradients, because the gas dynamics 
in the bulk depends, in particular, on the boundary conditions for the burnt matter. In the 
framework of the weak nonlinearity expansion, this non-locality shows itself as the necessity 
to increase the differential order of the equation for the flame front position. For example, 
it was shown in Ref. that in order to take into account influence of the vorticity drift on 
the structure of stationary flames with the accuracy of {6 — 1)^ in the asymptotic expansion 
for 6 —>■ 1 (which is the relevant weak nonlinearity expansion in the stationary case), one 
has to increase the differential order of the integro-differential equation by one as compared 
to the Sivashinsky equation (the latter is of the second order in 6' — 1, corresponding to 
the approximation in which the gas flow is potential on both sides of the front). Roughly 
speaking, the nonlocal relations, being integral from the non-perturbative point of view, are 
treated differential of inflnite order in the framework of the weak nonlinearity expansion. 

Let us now turn back to our present purpose of derivation of equation for the front po- 
sition at the second order of nonlinearity. Remarkably, it turns out that this approximation 
is exceptional in that the above-mentioned nonlocal complications do not arise in this case. 
We will give now a simple illustration of this important fact. Let us consider a weakly curved 
flame propagating in ^-direction with respect to an initially uniform fuel, and denote x the 
transverse coordinates. If the flame front position is described by equation z = /(x, t), then 
the condition that front is only weakly curved implies that c}/(x, t)/c}x is small, and so is 
the gas velocity perturbation 5v. From the continuity equation 



divv = 



and the Euler equations 



'dt 



+ (vV)v = --VP, 



P 



one has 
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V^P = -p diVkdkVi , (1) 
where summation over repeated indices is understood. Equation (|l|) implies 

P{z, X, t) = —p J dz'dx'G{z, X, z' , x', t)diVk{z', x', t)dkVi{z', x', t) + il{z, x, t) , (2) 

where f2(z,x, t) is a local function of the flow variables, satisfying V^fi = 0, and 
(^(z, X, z', x', t) is the Green function of the Laplace operator, appropriate to the given 
boundary condition. Note that the latter, being a condition on the pressure jump across the 
flame, is imposed only at the flame front. Since the unperturbed velocity field is spatially 
uniform, one has from Eq. @ for the curvature induced pressure variation 6P 

6P{z,:x-,t) = —p J dz'dx'G{z,x, z' ,x.\t)di6vk{z' jX.' ,t)dk6vi{z' jX.' ,t) + 6il{z,x,t) . (3) 

Next, let us denote Sx the characteristic length of the front perturbation. Then the corre- 
sponding length in ^-direction 

df 



Srr ■ 



Of 




df 


dt 




9x 



In other words, integration over z' in Eq. (|^) is effectively carried out for z' e 
(/(x, t), /(x, t) + 6z). Furthermore, as we will see in the sections below, 

|5v| 

Therefore, if one is interested in evaluating the pressure variation at the front, one can 
rewrite Eq. (^) with the accuracy of the second order 

(5P(/(x,t),x,t) = -p J dx'G{x,x',t)d,6vk{f{x',t),x',t)dk6v,{f{x',t),x',t) 

+ 5fi(/(x,t),x,t), (4) 
G'(x,x',t) = / dz'G{f{x,t),x,z',x.',t) . 



As we noted above, "boundary condition" for the pressure field is imposed only at the flame 
front. Therefore, the Green function G is independent of any other conditions relevant to 
the flow of the burnt matter [e.g., boundary conditions on the tube walls, in the case of 
flame propagation in a tube). The same is true for the function Q, since the value of f2 at a 
given point depends only on the value of flow variables at the same point. We thus see that 
Eq. (^) is an integral relation between local functions of the flow variables, defined on the 
front surface, which is independent of the flow dynamics in the bulk. Furthermore, since the 
right hand side of Eq. (1) is of the second order, it is not difficult to show, using the linear 
decomposition of the flow field into potential and vortex components, that 2;-derivatives of 
the velocity can be expressed via its derivatives along the fiame front, bringing this equation 
into the transverse representation. This implies that at the second order, there exists a 
universal equation which describes the fiame front dynamics in terms of the front position 
alone. In practice, it is actually more convenient to work with the differential form of the 
flow equations, rather than integral. How their transverse representation can be derived will 
be shown in detail in Sec. [11 A[. On the basis of this result, the nonlinear non- stationary 



equation will be derived in Sec. |ITI| . Its analysis is carried out analytically using the method 
of pole decomposition in Sec. The obtained results are summarized in Sec. 0. 
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II. FLOW EQUATIONS AND CONSERVATION LAWS 



Let us consider a 2D flame propagating in the negative z-direction with the normal 
velocity U{ with respect to an initially uniform quiescent combustible gas mixture. Denoting 
X the transverse coordinate, we introduce dimensionless space and time variables {i],C,) = 
{x / R, z / R) , T = tUf/R, where i? is a characteristic gasdynamic length of the problem. Then 
the normalized fluid velocity v = (v^/UfyVz/Uf) = {w,u) and pressure U = {P — = 
—oo))/p^Uf obey the following equations in the bulk 

du dw ^ 

du ^ du ^ du 1 dU 

dr dC, drj Q d^ ' 

dw dw dw 1 dli , . 

dr dC, drj Q dr] ' 

where g is the fluid density scaled on the density p_ of the fuel. 

The above flow equations are complemented by the following conservation laws to be 
satisfied across the flame front 

df 

u+-u_- ^{w+ - w^) = {6- 1)N , (8) 
w^-w. + ^{u^ -u_)=e\n9 (^Dw^ + ^Du_ + j^D^^ , (9) 

ehiO (d^f 2^ d^f ^^d^f ldfdN\ 
N \ dr'^ drdrj ~ d'rf N drj drj J ' 

where 



D = — +\w^ + —-^]—, N 




N dr] j drj ' \ ~^ \dri 



'df 



and e is the small dimensionless ratio of the flame thickness to the characteristic length. 
Finally, the following so-called evolution equation 

df df 9\ne fdN d , , d^f\ 

completes the above system of hydrodynamic equations and conservation laws, determining 
dynamics of the flame front itself. 

Below, we will need the general solution of the flow equations upstream. Going 

over to the rest frame of reference of the initially plane flame front, this solution is readily 
obtained as follows. Since the flow is potential at ^ = — oo (where m = 1, w = 0), it is 
potential for every < f{r], t) in view of the Thomson theorem 0, thus 
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u 



l + u = l+ J dk Ukexp{\k\^ + ikr]) , (12) 



w = Hu, (13) 
du $ 

where the Hilbert operator H is defined by 

(^/)(r/) = -p.v. / dC 

TX J 



^ + 4n + -(M= + u.=)=0, (14) 



i„,_ V" /(c) 



"p.v." denoting the principal value. Equations (|T2p, (|T3p represent the general form of 
the potential velocity filed satisfying the boundary conditions at ^ = — oo, while Eq. ( p!4[) 
is nothing but the Bernoulli equation. Note that the 2D Landau-Darrieus operator $ is 
simply expressed through the Hilbert operator $ = —OH. Note also that although the 
relation w = Hit between the velocity components upstream is nonlocal, it is expressed in 
terms of the transverse coordinate r] only. 



A. Bulk dynamics in transverse representation 

Our next step is the reduction of the system of fiow equations (^ - (|^) to one equation 
in which the role of the coordinate ^ is purely parametric. For this purpose, it is convenient 
to introduce the stream function ip via 

dip dip 

The stream function satisfies the following equation 

vV)V^^ = 0. (16) 




1. First order approximation 

To perform the weak nonlinearity expansion, it is convenient to explicitly extract zero- 
order values of the fiow variables downstream 

u = e + u, n = -e + i + n. (17) 

Then in the linear approximation, Eq. (|16D takes the form 
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(18) 



Its general solution can be written as a superposition of the potential and vorticity modes 
satisfying, respectively, 



General solution of Eq. (|T9D has the form analogous to Eqs. (|T2|), (|l^) 



+00 



(19) 
(20) 



(21) 



-Hup . 



Differentiating Eq. 



with respect to rj we obtain 
9 Or di] 



(22) 



(23) 



Next, linearizing Eq. (|^) and using the above equations for the potential and vorticity modes, 
one finds 



n, = 0. 

With the help of Eqs. (|3), (|2^) , and (|25|) equation (|^) can be rewritten as 

dw 1 du ^ 
dr] 9 Or 



(24) 
(25) 



(26) 



In this form, the flow equation governing dynamics downstream contains no explicit opera- 
tion with the ^-dependence of the flow variables. In other words, this dependence is rendered 
purely parametric. Let us now show that Eq. (|26|) can be generalized to take into account 
interaction of the perturbations. 



2. Second order approximation 
At the second order, Eq. (|16D takes the form 



d_ 



+ w 



(1) 



dr] 



(27) 
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General solution of this inhomogeneous equation is the sum of general solution of the ho- 
mogeneous equation, given by Eqs. (^)-(|25|), and of a particular solution ipj^^ which can be 
chosen to satisfy the following equation 



Written in components, Eq. ( pq ) has the form 

_^ + ,_^ + .<.)^_-_j =0, (29) 

Next, retaining the second order terms in the Euler equation @ and using Eqs. (|^), (p3|), 
(0) one can obtain the following relation between the velocity and pressure fields 

Finally, taking into account explicit structure of the potential mode and Eqs 
(^) it is not difficult to verify that Eq. (^) can be rewritten in terms of the sum 
+ + vi^^ as follows 



[w,u) 



du Bw i (~ u^ + w'^\ ( du Idw dll\ „ 
_-»_-*^OT + ^j+..^- + -- + -]=0. (32) 

which is the transverse representation of the flow equations at the second order of nonlin- 
earity we are looking for. 



III. NONLINEAR EQUATION FOR THE FLAME FRONT 

Now we can show that the set of equations (|]) - (PD , ([l3|) , (0), and (|32D can be reduced 
to one equation for the function f{ri,T). To this end, it remains only to rewrite the right 
hand sides of the conservation laws and evolution equation in the form in which transition 
on the flame surface — > f{ri,T)) is performed after all differentiations, again bringing 
the latter to the transverse form. This is easily done using the continuity equation and 
taking into account potentiality of the flow upstream 
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After having done this, one sees that the knowledge of explicit ^-dependence of the 
flow variables turns out to be unnecessary. Roughly speaking, the ^-dependence of a func- 
tion F{^,r],T) describing the flame front shape is known in advance, since the equation 
F{C:V:'^) = can always be brought into the form ^ — fiji^r) = (with / many-valued, 
in general). To determine the flame front evolution, therefore, it is sufficient to find only 
[r], r)-dependence of the functions involved. To put this intuitive reasoning in the formal 
way, it is convenient to introduce separate designations for the up- and downstream velocity 
and pressure fields, distinguishing them by the 5U]9erecripts " — " and " -|- " , respectively. 
Then, setting ^ = f{r],T), equations (|T3p, (0), ( [5^ ) together with the conservation laws (^ 
- (p!0|) and evolution equation (p!T|) can be rewritten identically as follows 



du 
dT 



W 

+ |,n- + f{(M-)2 + (^-)2} 



1 dw+ I afi+ \ 

e dr ~^ dri J 



Of] ' 



— W 



w 



u 



n+ -n- 



V 



u 



or] 



dl 
dr 



Hu 


6>-l 

.(0-l)0 + 5ln^0 



dlY 
9vJ 



9r,2 ) I 



(*) 



«=/{'?,t) 



Suppose we found a solution / = f{r],T), v~ = v' r), v+ = v+(,^, r), etc. of 
the set of equations in the large brackets in (*). Then, in particular, these equations are 
satisfied for ^ = f{ri,T). On the other hand, since no operation involving ^ appears in these 
equations, the function f{r],T) is independent of a particular form of .^-dependence of the 
fiow variables. For the purpose of deriving an equation for f{ri,T), it is most convenient 
to work with the particular solution in which all the functions are ^-independent, and to 
omit the large brackets in (*). Therefore, we can replace the above set of equations by the 
following 

(36) 
(37) 

(38) 



~~d7 



1 duj^ dn 




dr] 



U!~^ — to + Tr~(ti — f ) 
or] 



TT 



V 



dr] 



UJ 



■ TT 

dr 



l)^ + eln^^ 

e\ne (duj- 




+ 



9 — 1 \ dr] dr] 



(39) 
(40) 
(41) 

(42) 



where v,uj, and vr are the ^-independent counterparts of the fiow variables u,w, and 11, 
respectively, and 
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v~ = l + v~, v^ = e + v+. (43) 

The fact that now the function /(?], r) does not enter the arguments of these variables 
allows us to avoid expanding them in powers of /, employed in Ref. 0. In fact, such an 
expansion is irrelevant to the issue, since all the equations governing flame propagation are 
invariant with respect to the space translations, and therefore, all terms containing powers 
of undifferentiated / should appear in invariant combinations in the final equation for /. In 
view of this invariance, the function / itself does not need to be small even if the front is 
only weakly curved. We thus see that the /-dependence of the flow variables through their 
arguments must eventually cancel out in some way in any case. 

Now, the system of Eqs. ( PU]) - (^) can be transformed into one equation for the function 
f{ri,T). To simplify the derivation, the stabilizing effects due to the finite flame thickness 
will be taken into account only in the linear approximation. As we already mentioned in 
the Introduction, these effects are mainly taken into account in order to provide the short 
wavelength cutoff for the spectrum of the flame front perturbations. In view of this, the 
nonlinear corrections in the e-terms are of little interest. Using Eq. ([36|), the "evolution 
equation" (|42D can be rewritten within the accuracy of the second order 





Next, solving Eqs. (^) - (|4l| ) with respect to v~^, uj~^, vr"*", and substituting the results 
into Eq. (|38|) , one obtains an equation for the upstream variables v~ , uj~ , tt~, which can 
be further reduced to an equation for v~ alone using Eqs. (|36|) , (W^). Upon substituting 
the expression (|^) into the latter equation, one arrives at the nonlinear equation for the 
function /(r^, r), which we do not write explicitly because of its great complexity. It can be 
highly simplified using the first order LD-equation 

1^.4/. .^^(/^TTTi-l). (45) 

in the second order terms and in the terms containing e. The above value for u corresponds 
to the exponentially growing solution in the LD-theory. The nonlinear equation thus takes 
the form 

where 

(z/ + l)2 3z/ + l + 
a = — 77:^ \- 



4^2 4 40 

0-1 2 

^ AO 

6 = - , ^ ly + e] InO 

e-1 \ 2 2 2 
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Following Zhdanov and Trubnikov Eq. ( pq ) can be further simplified by rewriting its 
linear part in the form 

Since we are only interested in the development of unstable modes of the front perturbations, 
which satisfy Eq. (^Sf), we can transform the first factor in the right hand side of Eq. 
as follows 

6 + 1^ , d 6-1 6 + 1 6-1 

\ > H = 

26 dr 2z/ 26 2u 

Therefore, Eq. (^) becomes 

f--/--'J^g)^-KI)^^^(*^^^^i^»- '''' 

where 

a,/3,7,(5 = , 

Finally, we would like to comment on the range of validity of the derived equation. Gen- 




erally speaking, Eq. (48) is only applicable for description of the early stages of development 
of the LD-instability, since it is obtained under the assumption of smallness of the front 
slope. Even if the fiame evolution is such that it smoothly ends up with the formation of 
a stationary configuration (instead of spontaneous turbulization), this assumption becomes 
generally invalid whenever the process of fiame propagation is close to the stationary regime. 
In fact, it can be easily shown that the assumptions of stationarity and weak nonlinearity 
contradict each other (detailed discussion of this point can be found in Ref. ||^). Inciden- 
tally, the fact that the transition to the stationary regime in Eq. (|48|) is formally incorrect 
is clearly seen from its derivation given above. Namely, the stationary form of this equation 
depends on the way the first order relation (^) is used in the second order terms before 
time derivatives are omitted. Only in the case of small gas expansion {6 — 1) is weak non- 
linearity approximation justified at all stages of development of the LD-instability, in which 
case Eq. (^) goes over to the well-known Sivashinsky equation 

9^ + 2(9^) ^^^^'^W' ^^^^ 

since 

P = l + 0i{6-iy), 7 = O((0-lf), 5 = -l + Oi6-l), u=^-^ + 0{{6-iy). 

In this respect, a natural question arises as to what extent equation (|^) is actually valid 
when 6 is arbitrary. Since the structure of higher order terms of the power expansion is 
unknown, it is very difficult to give even a rough estimate. Leaving this question aside, we 
will simply assume in what follows, that this equation is formally valid for all times. It will 
be shown below that at least in the case of fiame propagation in narrow tubes, solutions to 
the stationary version of Eq. (^8]) are in reasonable agreement with the results of numerical 
experiments [0] for fiames with the gas expansion coefficient up to 6* ~ 3. 
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IV. THE POLE DECOMPOSITION 



As in the case of ZT-equation, development of the LD-instabihty of a plane flame can 
be described in terms of the pole dynamics. To show this, we first perform the following 
nonlinear transformation 



/ 



In terms of the new function 0, Eq. 
order. 



(50) 



takes the form, within the accuracy of the second 



dr 



v^cj) +{13 + va) 



+ 7 



e6 



0. 



(51) 



Spatially periodic (with period 26) solutions of equations of the type (|5T|) can be found 
using the following pole decomposition P,pi? 



2P 



(p{ri, t) = 0o(i") + a ^ In sin 
fc=i 



TT 

26 



(52) 



where the value of the amplitude a as well as dynamics of the complex poles rik{T), k = 
1,...,2P are to be determined substituting this anzats into Eq. (pT|). Since the function 
0(?7, r) is real for real rj, the poles come in conjugate pairs; P is the number of the pole 
pairs. Requiring the 26-periodic solutions to be symmetric with respect to the reflection 
1] —7], one can obtain periodic as well as non-periodic solutions to Eq. ( p]) in the domain 
rj G (0, 6), satisfying the conditions |^(0, r) = f^(6, r) = 0, describing flame propagation in 
a tube of width 6 with ideal walls. 
Using the formulae H 



dr] 

cot X cot y 



vra 
26 



^ < 1 + 2 sign(Im rjk) cot 
k=i ^ 

- 1 + cot {x — y) (cot y — cot x) 



26 



sign(xj 



X 

\x\ 



(53) 



it is not difficult to verify that Eq. (|5T|) is satisfied by (f){ri, r) taken in the form of Eq. 
provided that 



and the poles rikir), k 
dr]k 



a = eSx, 
d(j)o 



dr 

, 2P, satisfy 



xicr'P' - uaP) , 



dr 



+ i sign(Im rjk) (z^ + 2P^ax) 



2P 



-^X {7sign(Im rj^ Im rjm) - (3 - va} cot 



m=l 



TT 

26 



iVk - Vn 



(54) 
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where the following notation is introduced 



a 



>0, x={.l3 + T^a-iy 



Since the application of pole decomposition to Eq. (|5l|) is quite similar to that given in 
Refs. we will present below only final results, referring the reader to these works for 



more detail. 

Following Ref. [|I^, we first consider two poles {r]i,ri2) in the same half plane of the 
complex 7], which are fairly close to each other, so that their dynamics is unaffected by the 
rest. Then one has from Eq. (^) 

d , , Ae6 



dr ?7i - ?72 ' 

which indicates that the poles attract each other in the horizontal direction (parallel to the 
real axis), and repel in the vertical direction, tending to form alignments parallel to the 
imaginary axis. Furthermore, assuming that the pole dynamics ends up with the formation 
of such a "coalescent" stationary configuration, and using the fact that 7 < /5 + z/a (it is 
not difficult to verify that actually 7/(/3 + z/a) < 1/3), the following upper bound on the 
number of pole pairs P can be easily obtained from Eq. (|54D 

- 2 2a 

Still, for sufficiently wide tubes (such that a < i^/S), the solution (|52D is not unique: different 
solutions corresponding to different numbers P of poles are possible. To find the physical 
ones, it is necessary to perform the stability analysis. Noting that the functional structure of 
Eq. (^) is very similar to that of Eq. (^9[) , the stability analysis of Refs. |ll|]- |jl3|) where it 



was carried out for the Sivashinsky equation, will be carried over the present case. According 
to this analysis, for a given not-too-wide tube, there is only one (neutrally) stable solution. 
This solution corresponds to the number of poles that provides maximal flame velocity, i.e., 



Pr^ = Int 



u 1 
2^^ 2 



Int(x) denoting the integer part of x. Thus, the flame velocity increase Wg = {dcpo/drl of 
the stable solution can be written as 

where 

2 

W^ = ^. (56) 

It is seen from Eq. ( ^T]) that the spectrum of front perturbations is effectively cut off at 
the wavelength 

A = H=i^. (57) 
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representing the characteristic dimension of the flame cellular structure. 

Fig. ID compares the theoretical dependence of the maximal flame velocity increase 
on the gas expansion coefficient, given by Eq. (|56D, with the results of numerical experiments 
1^ . For comparison, we show also the corresponding dependence calculated on the basis of 



the Sivashinsky equation. Dependence of the effective cut-off wavelength on the expansion 
coefficient is shown in Fig. ^ We see that even beyond of its range of applicability, Eq. (|5TD 
provides reasonable qualitative description of ffames with the expansion coefficients of prac- 
tical interest. Complete investigation of the LD- instability on the basis of this equation will 
be given elsewhere. 



V. CONCLUSIONS 

The main result of our work is the nonlinear non- stationary equation ( ^6|) which describes 
development of the Landau-Darrieus instability in the second order of nonlinearity. We have 
derived this equation on the basis of the only assumption of smallness of the flame front slope. 
Thus, nonlinear evolution of the front perturbations generally obeys Eq. which takes 
even simpler form (|5l|) if one is only interested in dynamics of the exponentially growing 
LD-solutions. It is important to stress that no assumption concerning the value of the flame 
front position has been used in the derivation. Therefore, Eq. (^) can be applied not only 
to plane flames, but also to the problem of unstable evolution of any flame conflguration, 
provided that the front slope is small. 

It is also worth of emphasis that Eq. (^) is obtained without any assumptions about the 
structure of the gas flow downstream. Thus, this equation is the direct consequence of the 
exact hydrodynamic equations for the flow flelds in the bulk, and conservation laws at the 
flame front. We would like to stress also once again that the universal form of Eq. ( ^6] ) is 



the distinguishing property of the second order approximation. In the general case, equation 
for the flame front should contain also information about the flow of the burnt matter in 
the bulk. Indeed, as was shown in Ref. 0, the boundary conditions for the burnt matter 
are invoked in the course of derivation of the equation already at the third order. This 
universality of Eq. P6| ) allows it to be widely applied to the study of flames with arbitrary 
front conflguration, propagating in tubes with complex geometries. 
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FIGURES 




FIG. 1. Maximal flame velocity increase Wm versus the gas expansion coefficient, given by 
Eq. (|56[); the marks are according to Ref. [8]. Accuracy of the experimental results is about 20%. 
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FIG. 2. The cut-off wavelength scaled on the flame thickness versus the expansion coefficient, 
^iven by Eq. (p^). 
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